Exchange Coupling in a One-Dimensional Wigner Crystal 
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We consider a long quantum wire at low electron densities. In this strong interaction regime a 
Wigner crystal may form, in which electrons comprise an antiferromagnetic Heisenberg spin chain. 
The coupling constant J is exponentially small, as it originates from tunneling of two neighboring 
electrons through the segregating potential barrier. We study this exponential dependence, properly 
accounting for the many-body effects and the finite width of the wire. 
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Quantum wires exhibit a plethora of interesting phe- 
nomena. In particular, quantization of conductance, 
which is a fundamental manifestation of the quantum na- 
ture of the electron, has been actively studied ever since 
its first observation in quantum point contacts Q. The 
phenomenon presents itself as very flat plateaus of linear 
conductance G at integer multiples of Go = 2e 2 /h, as 
a function of gate voltage which tunes the electron den- 
sity in the wire. Since the first observation of the phe- 
nomenon, its various facets have been studied by mea- 
surements of thermal transport 0, Q , noise , and ex- 
periments on systems with superconducting elements . 

A new generation of experiments in quantum wires has 
revealed an unexpected structure at low electron densi- 
ties: a plateau at about 0.7 Go for short E i pl [ § , I ^ Hf H and 
at 0.5 Go for long quantum wires [TlL fT3 fl3L Il4| . These 
new features have generated much interest as they are 
likely caused by electron correlation effects. The origin 
of the new plateau has not yet been established; how- 
ever, the experiments 0, 0, 13 point to the important 
role played by the electron spins. 

In recent experiments of a different kind, involving tun- 
neling between two parallel ballistic quantum wires, map- 
ping of the spectrum of spin and charge excitations was 
achieved dll- In addition, unexpected behavior indicat- 
ing electron localization was observed at low densities. 
Of particular interest is the concurrence of the localiza- 
tion with the drops in the conductance steps, which in- 
dicates Il5la possible connection with the 0.7 structure 
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Recent theoretical work suggests that qualitatively 
new transport properties of quantum wires at low elec- 
tron densities may be due to the formation of a Wigner 
crystal state of electrons. The latter is expected to occur 
when the density is low enough for the potential energy 
of the Coulomb repulsion to overwhelm the kinetic en- 
ergy of the electrons in the wire. In the crystalline con- 
figuration, the electron spins form an antiferromagnetic 
Heisenberg spin chain, with the exchange coupling con- 
stant J emerging as a new energy scale. The Heisenberg 
exchange can be viewed as arising from tunneling of two 
neighboring electrons through the potential barrier cre- 
ated by their mutual repulsion, and by the repulsion from 



all other electrons of the wire. As a result, J is expected 
to be small compared to the Fermi energy Ep, and qual- 
itatively new transport properties of the quantum wires, 
such as the plateau at 0.5 Go, are expected 0, 0] at 
temperatures in the range J <C T <C E p. In addition 
to conductance measurements |ll|, |l2j, |l3|, Uj£ , the small- 
ness of J should have a strong effect on the velocity of 
the spin excitations in the wire measured in Ref. and 
on transport properties of quantum wires in the in-plane 
magnetic field. Therefore it is important to have a reli- 
able estimate of J. 

In this article we evaluate the exchange coupling con- 
stant J in a quantum wire at low electron density, when 
the electrons are in a Wigner crystal state. Similarly to 
the case of a two-dimensional crystal [lijj . the Coulomb 
interaction couples other electrons of the wire to the tun- 
neling of the exchanging pair, thus turning it into a many- 
body process. We study this process, using a combina- 
tion of analytic and numerical tools in two distinct ap- 
proximations. First we consider a wire of zero width, and 
subsequently account for the effect of a finite width. 

Disregarding the collective nature of the tunneling pro- 
cess, an estimate of J can be made by considering 
a single pair of exchanging electrons, with all the others 
held in fixed positions separated by distance b = n . 
The only dynamical variable in this case is the relative 
coordinate x of the two electrons. Tunneling through the 
potential barrier U(x), created by their mutual repulsion 
and by the repulsion from all the other electrons, lifts the 
ground state degeneracy caused by the symmetry with 
respect to inversion x — > —x. The exchange energy J 
coincides with the exponentially small level splitting. In 
the semiclassical approximation it is given by 



J = J* exp ( — -^j , So = J Axy/ mU(x), 



(1) 



where m is the effective mass of the electrons. 

Determining J for a given density amounts to calculat- 
ing the exponent Sq/H, with the prefactor J* being of sec- 
ondary importance for the problem at hand. The action 
in Eq. Q is easily brought to the form So — firjo/y/naB, 
and the constant 770 ~ 2.817 is evaluated numcr ically [Hi. 
(Here e is the dielectric constant, and as — eh 2 /me 2 is 
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the effective Bohr radius of the semiconductor host) . The 
precision of this estimate is unclear, as only the two ex- 
changing electrons were allowed to move. 

To account for the motion of all the electrons, we 
treat the tunneling using the instanton method, in which 
the action is represented as an imaginary-time inte- 
gral. By measuring distance and time in units of b and 
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respectively, we find the dimensionless ac- 
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Xj(T)-Xi{T) 



(2) 

Here Xj is the dimensionless coordinate of the j-th elec- 
tron. In this approach the exchange constant is given by 



J = J* exp 



(3) 



where r\ is the action (0) minimized over the trajectories 
Xj(r) of all the particles. 

The approximation Q is equivalent to imposing a con- 
straint Xj(r) = j for all j ^ 0, 1. By releasing this con- 
straint we minimize the action over more variables. Thus 
allowing for the participation of all the electrons in the 
exchange process will result in r\ < t]q. 

In the course of the tunneling, near neighbors of the ex- 
changing pair undergo displacements comparable to the 
inter-electron distance. The contribution of these elec- 
trons to the tunneling action So can only be studied by 
numerical means. By contrast, the displacements of re- 
mote electrons are small and vary smoothly, allowing for 
a continuous description and analytic treatment of their 
contribution to r\. Combining the two contributions per- 
mits a reliable evaluation of r\. 

We study the long distance contribution first. To this 
end, we modify the problem as follows. Starting with the 
nearest neighbors of the exchanging pair, a large number 
I ^> 1 of electrons are assumed to be constrained to their 
equilibrium positions. As a result, all the moving neigh- 
bors are far from the exchanging pair, which enables an- 
alytic treatment of the problem. Furthermore, it will be 
convenient to assume that only TV 3> I of the subsequent 
electrons can move. Apart from taking the limit N — > oo, 
this assumption will enable us to study how this limit is 
approached at large N. 

We treat the N mobile electrons as a continuous string 
described by its displacement u(X,t) from equilibrium, 
I < X < I + N . Since I ^> 1, the string is far from the 
exchanging pair, and the resulting displacement is small, 
u <C 1. This enables us to account for the coupling of 
the string to the exchanging pair in linear order in u and 



present the string contribution rj s to the action as 



duj x 
2^^ 



X 4 



]k denotes the expansion coefficient of 
4 in the eigenmode basis. The coupling 



(4) 

Here the parameter k labels the (normalized) cigenmodes 
\[ / fc(X) of the string, Uk is the respective plasmon fre- 
quency, and [X 
the function X' 
of the string to the exchanging pair is accounted for 
by the function f u defined as the Fourier transform of 
/(r)=JC (r)X 1 (r). 

According to Eq. (@J the coupling of a remote mobile 
electron to the exchanging pair decays as 1/A 4 . This 
can be understood by noticing that symmetric displace- 
ment of the exchanging electrons, X\ — 1 — Xq, creates a 
quadrupole potential cx 1/A 3 at large distances. To lin- 
ear order in u(X), the displacement of a mobile electron 
is equivalent to creation of a dipole at point X, whose in- 
teraction with the quadrupole field decays as 1/A 4 . Since 
X > I, the coupling of the string to the exchanging pair 
is weak at I ^> 1 . The feedback effect of the string upon 
the exchanging pair is further reduced by the same pa- 
rameter. Thus we neglect the latter effect and substitute 
into Eq. the function f u computed with unperturbed 
functions X (t) and Xi(t). 

We then minimize the action (0J, and find the string 
contribution to the total action 
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Here we used the fact that at frequencies well below the 
Debye frequency of plasmons the Fourier transform f u of 
fir) smoothly approaches its zero- frequency limit fo- 

The eigenmodes ^k(X) are essentially standing plas- 
mon waves in the string with wave number k; only near 
the edges does their shape slightly deviate from the sinu- 
soidal one. Using the proper plasmon dispersion 



= fcV21nfc-i 



(6) 



we find the following asymptotic behavior in the limit of 
large N: 



Vn =V + 



a 



(7) 



The expressions for the two constants r\ and a involve 
parameters fo and ^' k (l) which are easily found numeri- 
cally. As a result we obtain 



v = m 



0.020 
K/InT 



0.011 



(8) 



In the limit of large number of mobile particles N — > 
oo, the numerator of the exponent in the expression @ 
for the exchange constant takes the limiting value r\\ as 
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FIG. 1: Triangles and circles represent the calculated expo- 
nents. Only the values represented by the circles (N — 6 — > 
11) were used to fit to the analytically predicted asymptotic 
behavior 0. The horizontal line indicates the value of the 
exponent in the limit when all electrons participate in the 
tunneling. For this particular fit we obtained r\ = 2.7981. 



expected, 77 < r/Q. Somewhat surprisingly, however, the 
total effect is numerically small: even at I ~ 1 we find 
Vo ~ V ~ 0.02, i.e., the total correction is of the order of 
one percent of the unperturbed value 770 ps 2.817. 

It is important to note that the above approach re- 
lies on the large parameter / 3> 1, which is the number 
of immobile electrons between the mobile ones and the 
exchanging pair. In the actual exchange process, all elec- 
trons move and contribute to the action. Fortunately, 
the contribution to the tunneling action from a few near 
neighbors of the exchanging pair can be found numer- 
ically. To this end, we minimize the action (0) with 
a finite number N of mobile electrons on each side of 
the exchanging pair, while keeping the remaining elec- 
trons at their equilibrium positions. Mathematically, 
this procedure is equivalent to solving a system of N 
second order differential equations of motion with the 
boundary conditions Xi(— 00) = 1, Xi(+oo) — 0, and 
Xj(-oo) = Xj(+oo) = j for j = 2, 3, . . . , N + 1. In 
practice, we were able to perform this calculation for up 
to N = 10; the results for the respective action t)n are 
shown in Fig.^ The calculated asymptotic behavior Q 
can be used to fit these results. For instance, using the 
last six calculated values 775, ... , 7710, we find an excellent 
fit with fitting parameters 77 = 2.7981, and a — 0.0221. 
In order to obtain a more accurate result and estimate 
uncertainty bounds, we considered a series of two differ- 
ent kinds of fitting schemes. One consisted of fits using 
various numbers of 77^ points; for the other, we intro- 
duced a third fitting parameter, namely an offset to N. 
Note that the latter does not contradict the predicted 
asymptotic behavior 0. By combining the analytical 
treatment with numerical calculation in this fashion we 



have obtained the exponent 77 = 2.79805 ± 0.00005 in 
Eq. @ in the limit when all the electrons participate in 
the exchange process. 

Since the effect of neighboring electrons upon 77 in 
Eq. @ is fairly weak, it is important to consider other 
corrections which may potentially affect the exchange 
coupling in quantum wires. The most important cor- 
rection is due to the fact that quantum wires have finite 
width w. As a result, the exchanging electrons can utilize 
the dimensions transverse to the wire axis to minimize 
the Coulomb repulsion at short separations. Therefore, 
the short distance behavior of the interaction is strongly 
affected by the regularization of the singularity inherent 
in the previous one-dimensional idealization. The impli- 
cation is that the action should be further reduced. 

There are two limits for which we can calculate the 
finite-width corrections to the action. The small width 
limit, w <C as, implies that the typical electron inter- 
action energy e 2 few is much smaller than the subband 
spacing in the wire H 2 /mw 2 . Then the processes of elec- 
tron scattering to higher subbands are negligible, and 
the effect of the finite width of the wire amounts to 
small smearing of electron density in transverse direc- 
tion. Thus we calculate an effective one-dimensional in- 
teraction U w (x) as an average of the three-dimensional 
Coulomb interaction over the transverse coordinates, us- 
ing the ground state wave function in the confining po- 
tential. This procedure changes the interaction potential 
in Eq. JIJ only at distances \x\ < w, which enables us to 
estimate the correction to the action Sq as 



AS = I dx 

-b 



[yj mU w (x) — \J mU(x)j = —hn^ 



(9) 

where only the numerical prefactor k is sensitive to the 
details of the confining geometry. As an illustration, we 
consider an axially symmetric quadratic confining poten- 
tial and define w as root-mean-square deviation of the 
electron from the axis in the ground state. At short dis- 
tances I a; I <C b, we find an effective interaction 



U w (x) = 



2we 



1 



- 1 £ 



(10) 



which, of course, coincides with e 2 /e\x\ at \x\ 3> w. Sub- 
stituting this expression into Eq. ©, we find k ~ 2.99. 
We estimated the corrections to AS from small admix- 
ture of higher subbands to be of order (w/as) 5 ^ 2 . 

In the large width limit, w 3> as, we can once again 
consider a semiclassical approximation analogous to the 
one-dimensional one we employed previously. Indeed, in 
the instanton approach the trajectories of two exchang- 
ing electrons behave as shown in Fig. The electrons 
pass each other at a distance ro, at which the Coulomb 
repulsion is of the order of the confining potential. For 
example, in the case of quadratic potential we estimate 
e 2 /ero ~ (h 2 /mw 2 )(ro/w) 2 and find r ~ (w 4 / 'cle) 1 ^ ■ 
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FIG. 2: The classical exchange trajectories obtained by solv- 
ing numerically the equations of motion for two electrons (zero 
density limit) in a wire of width w S> as with a quadratic con- 
fining potential. Arrows indicate the direction of motion. The 
unit of length is ro = 2a_B(w/as) 4 ^ 3 - 
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classical equations of motion in the inverted potential. In 
the low-density regime, n <C r^" , the trajectories deviate 
from straight lines only at short distances ~ ro, Fig. 
We thus evaluate the correction to So, by finding the 
difference of the two corresponding actions calculated up 
to a large distance cutoff. This calculation gives k 2.03. 

We are now in a position to calculate the exchange 
constant as a function of density for both the one- 
dimensional wire, and wires of various representative 
widths. The two corrections to the action, Eqs. © and 
have comparable magnitudes at w ~ Was- Consid- 
ering that for typical confinement geometries in realistic 
devices w does not exceed as, we use the finite- width 
correction in the form Jj|J. In Fig. we show the corre- 
sponding curves as a function of density. 

Summarizing, we have found the exchange energy in a 
Wigner crystal of electrons in a quantum wire in the ex- 
ponential approximation. The tunneling process involves 
many electrons neighboring the exchanging pair. Their 
effect was accounted for by combining a numerical calcu- 
lation for the near neighbors of the exchanging pair with 
an analytic calculation of the asymptotic behavior for 
the distant ones. Finally, we found the substantial cor- 
rections due to a finite width of the wire in both small- 
and large-width limits. 

Work supported by the U. S. Department of Energy, 
Office of Science, under Contract No. W-31-109-ENG-38. 
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FIG. 3: The calculated action as a function of density for var- 
ious wire widths w. Inset: the exchange energy J in Kelvins 
as a function of density for the same set of wire widths as in 
the main figure, using the estimate of the prefactor obtained 
in Ref. Il7l and parameters pertaining to GaAs. 



Since r Q w, the quantum particles always stay close to 
the instanton trajectory, and the semiclassical approach 
is applicable. Furthermore, since the two electrons pass 
each other at distance of order ro, their interaction po- 
tential is regularized at that length scale. This reduces 
the tunneling action in a way similar to the small width 
case ||SJ|, except that instead of w one should substitute 
r . Thus the correction to the exchange action is 



AS 



-hk ( — 

a B 



2/3 



(11) 



at ra > as. The correction grows faster as a function of 
w than in the regime w <C as- This generic feature of the 
problem occurs for any shape of the confining potential 
(other than the hard- wall one). However, the 2/3 power- 
law dependence is valid only for quadratic potential. 

To find the coefficient k, we compute the instanton for 
the exchange of two particles numerically by solving the 
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